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ABSTRACT 

We present the first metallicity gradient measurement for a grand-design face-on 
spiral galaxy at z ~ 1.5. This galaxy has been magnified by a factor of 22 x by a 
massive, X-ray luminous galaxy cluster MACS Jl 149.5-f2223 atz = 0.544. Using the 
Laser Guide Star Adaptive Optics aided integral field spectrograph OSIRIS on KECK 
II, we target the Ha emission and achieve a spatial resolution of Of'. I, corresponding 
to a source plane resolution of 170 pc. The galaxy has well-developed spiral arms 
and the nebular emission line dynamics clearly indicate a rotationally supported disk 
with Vrot/cr ~ 4. The best-fit disk velocity field model yields a maximum rotation 
of Vrotsini = 150 ± 15 km s^^, and a dynamical mass of M^yn = 1-3 ± 0.2 x 
10^° csc^(^) Mq (within 2.5 kpc), where the inclination angle z = 45 ± 10°. Based on 
the [Nil] and Ha ratios, we measured the radial chemical abundance gradient from 
the inner hundreds of parsecs out to ~5 kpc. The slope of the gradient is —0.16 ± 0.02 
dex kpc~^, significantly steeper than the gradient of late-type or early-type galaxies in 
the local universe. If representative of disk galaxies at 2; ~ 1.5, our results support an 
"inside-out" disk formation scenario in which early inf all/collapse in the galaxy center 
builds a chemically enriched nucleus, followed by slow enrichment of the disk over 
the next 9 Gyr. 

Subject headings: galaxies: abundances — galaxies: evolution — galaxies: high-redshift — 
gravitational lensing: strong 



1. Introduction 

The spatial distribution of heavy element abundances is a unique tool t o trace the stru cture 



formation history of disk galaxies. Since the very early work of ISearlel (| 1 97 1|) . IShieldsl (| 1 974|) . and 
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Zaritsky et al.l(|1994|) . it has been well established that there is a negative radial metallicity gradient 
in spiral galaxies such that the central disk region is more metal-enriched than the outer regions. 
The evolution of this gradient over cosmic time provides a powerful constraint on disk formation 
and evolution models. 

Sophisticated chemical evolution models(CEMs) have been successful in reproducing the cur- 
rent gradient of local gal axies, but different models disagree on the time variation of t he radial 



20091 : 



Marcon-Uchida et a 



time (Ichiappini et al.l 



1997 



Some models predict that the radial gradients steepen with 



metal l icity gradient (e.g.. iMoUa et all Il996l: IChiappini et al.l 1 19971 : iMagrini et all l2007l : Fu et al. 



200 lb . while oth er models predict that the rad ial gradients flatten with 



time (MoUa et al.ll 19971 : IHou et al.ll2000 : Prantzos & Boissiei]l2000l : lFu et al. 2009) . The reason for 
these opposite predictions is that CEM models are very sensitive to the prescriptions of the detailed 
physical processes that lead to the chemical enrichment of inner and outer disks. Observational 
constraints on these physical processes are lacking. 

Metallicity gradient evolution is unknown due to the small angular sizes and the strong decline 
of surface brightness of distant galaxies. Before the next generation of telescopes, the only way to 
overcome this observational hurdle is to combine the magnification power of gravitational lensing 
with the Laser Guide Star Adaptive Optics (LGSAO) aided integral field spectrograph units (IFUs ) 
on the w orld's larges t telesc opes (|jones et al.l 1201 Obi: IStark et al.l l2008l : ISwinbank et al.l |2007|) . 
Recently, Ijones et al.l (1201 Oa|) reported the metallicity gr adient of a strongly lensed dispersion- 
dominated galaxy at high-z ("the Clone arc" at z = 2.00). Ijones et al.l (|2010al) measured the slope 
of the gradient for the inner ~ 1 kpc and found it to be considerably steeper than local disk galaxies. 
The steep gradient at high- ^ is consistent with the "inside-out" disk formation scenarios. However, 
in a low resolution study, ICresci et al.l (1201 Ol) reported "inverted" metallicity gradients for three 
non-lensed 2; ~ 3 Lyman-break galaxies (LBGs), arguing for a "cold-flow" of primordial gas to 
the galactic core regions. Clearly, a large sample of galaxies with chemical abundance gradients 
across cosmic time are required to place robust constraints on disk formation theories. 

In this Letter, we report the first metallicity gradient measurement of a lensed spiral galaxy at 
z = 1.49 based on our OH- Suppressing Infra-Red Imaging Spectrograph (OSIRIS) observations 
at KECK II. Throughout this Letter we use a ACDM cosmology with Hq= 70 km s^ Mpc^, 
fiAf=0-30, and Qa=OJO. At z = 1.5, 1 arcsec corresponds t o 8.5 kpc and a look -back time of 9.3 
Gyr. We use solar oxygen abundance 12 -1- log(O/H)0=8.66 (|Asplund et al.ll2005h . 
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2. Observations and Data Reduction 



The spiral galaxy is a four-image lensed system first identified by ISmith et alj (120091) . lo- 
cated behind the massive, X-ray luminous cluster MACS Jl 149.5+2223. We observed the largest 
of these four imag es at (02000, ^2000) = (ll''49'" 35*^.284, -i-22°23'45'.'86) (called Spll49 hereafter). 
Smith et al.l (|2009h construct a detailed mass model for the galaxy cluster and show that Spll49 is 
one of the largest recorded lensed images of a single galaxy at z > 1. Owing to the lensing magni- 
fication and fortuitous face-on orientation, this galaxy exhibits more than 10 locally resolved H II 
regions (Figured), ideal fo r spatial chemical abundance studies. The star formation rate from the 
^555 band is ~ 6 M© yr'^ dSmith et al.ll2009l) . 
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Fig. 1.— Left: HST F814W, F555W-band color image of Spll49 behind the lensing cluster 
MACS Jl 149.5-^2223. North is 126° counterclockwise to the y+ direction. Spll49 is the largest 
of the four multiple images of a spiral galaxy at z = 1.49. The white box shows the 4'.'8 x 6'.'4 FOV 
of OSIRIS Hn3 filter with the (X'l per pixel plate scale. Right: wavelength-collapsed image of the 
Ho; emission from our OSIRIS datacube. Contours in black denote observed Ha flux levels of 0.6, 
1.3 in units of 10~^^ ergs s^ cm^ sec^. Labels on the left and bottom axes indicate the imaging 
plane coordinates while the right and top axes show the corresponding physical scale on the source 
plane. The white circles mark the three annuli (rl, r2, r3) we use for deriving metallicity gradient 
(Section HI). The Ha morphology aligns generally well with the rest-frame UV morphology. 



We observed Spl 149 OSIRIS (ILarkin et al 



2006|) in conjunction with the LGSAO system on 
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KECK n telescope on 2010 March 3rd. The weather conditions were moderate with the optical 
FWHM ranging between 0'.'6 and I'.'S. The average tip/tilt corrected FWHM during the observa- 
tion was O'.'l. To capture Ha and [Nil], we used the narrow band Hn3 filter (spectral coverage: 
1.594— 1.676/im; spectral resolution: ~3400), and a plate scale of O'.'l per lenslet. The observations 
were conducted in an ABAA dithering sequence with a lO'.'V West-North chop onto an object-free 
sky region. A total of 19 exposures were obtained for the target, with 900 s of exposure for each 
individual frame. The net on-target exposure time is therefore 4.75 hr. Figured] shows the OSIRIS 
field of view (FOV) and locations of Spl 149 on an HST F814W, F555W color image. 



Individual exposures were first reduced usin g the OS I RIS d ata reduction pipeline (ILarkin et al. 



20061). We used the sky subtraction IDL code of iPaviesI (|2007|) . which we have modified to opti 



mize the sky subtraction specifically around the wavelength region of Ha and [N ll]. Since we are 
mainly concerned with the emission lines, a first order polynomial function was fit to each spatial 
sample pixel (spaxel) to improve the removal of lenslet to lenslet variations. The final datacube 
was constructed by aligning the sub-exposures with the centroid of the Ha images and combing 
using a 3a mean clip to reject cosmic rays and bad lenslets. Telluric correction and flux calibration 
was performed by observing a few AOV standard stars immediately after the science exposure. 
The typical uncertainty in flux calibration is ~10%. 

Gaussian profiles were fitted simultaneously to the three emission lines: [Nil] A6548, 6583 
and Ha. The line profile fitting was conducted using a minimization procedure which takes 
into account the greater noise level close to atmospheric OH emissions. The centroid and velocity 
width of [Nil] AA6548, 6583 lines were constrained by the velocity width of Ha A6563, an d the 
ratio of [N ll] A6548 and [N ll] A6583 is constrained to the theoretical value (I Osterbrocklll989l) . An 



instrumental profile of 5.9A calculated from sky lines has been subtracted from the line widths. 
We demand a minimum signal-to-noise (S/N) of 3 to detect each emission line and a minimum 
S/N of 5 for kinematic and metallicity analysis. 



3. Kinematics and Lensing Modeling 

To align the astrometry to the IFU datacube, we first assign the brightest pixel of the Ha 
narrow-band image with the coordinates of the galactic core from the Hubble Space Telescope 
(HST) image, and then rotate the IFU datacube with the observing position angle. The Ha image 
can then be aligned with the HST image by matching coordinates. The upper right panel of Fig- 
ure \T\ shows our observed Ha emission line image, collapsed in wavelength, and smoothed with 
a Gaussian kernel of FWHM = 2.5 pixels. The Ha morphology aligns generally well with the 
rest-frame NUV/t/-band image, especially some of the H ll regions in the west. The SFR and size 
relation of individual H ll regions will be reported in a following paper (R. C. Livermore et al. 
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201 1, in preparation). 



We fit the velocity field with rotation disk models as described in lJones et alj(l2010b|) . The best 
fit yields a rotation velocity of Kot sin i = 150 ± 15 km s^^, with inclination angle i = 45° ± 10°. 
Using Mdyn = V^lt R/G, the derived dynamical mass is M^yn = 2.6 ± 0.4 x lO^^M© within R=2.5 
kpc. Spl 149 is a rotation dominated disk with Vrot/ cr ~ 4. 

To reconstruct the source-plane morphology and obtain flux magnification for Spl 149, we 



apply the most recent lens model for cluster MACS Jl 149.5+2223 as described in ISmith et al 



(|2009|) . The derived luminosity weighted magnification for Spl 149 is /i = 22 ± 2. Accounting for 
the lensing amplification, the intrinsic magnitude in the /—band is / ~ 23.4 ± 0.3, corresponding 
to Mb ~ -20.7. 

Spl 149 is stretched ~ 5 on each side {x,y) of the two-dimentional image. Because of this 
fortuitous uniform magnification, we can work with the observed datacube and divide a factor of 
5 linearly and a factor of 22 in flux to recover the intrinsic linear geometry and luminosity of the 
galaxy. Our results are therefore unchanged when we use the reconstructed source-plane datacube. 



4. Metallicity Gradient 



For the metallicity analysis in this work, we use the empirical stro ng line diagnostic 
N2 = log([N II] A6583/Ha) index, as calibrated bv lPettini & Pagell hoo4. hereafter PP04): 



12 + log(0/H) = 8.90 + 0.57 x N2 



the 



(1) 



We note that the N2 ratio is sensitive to shock excitation, the ionization state of the ga s , and the hard 



ionization radiation fi eld from an active galactic nucleus (|Kewley & Dopitd l2002l : iKewley et al 



2006l : lRich et al.ll2010|) . Our N2 ratios are consistent with star-formation throughout the disk. The 
ratio of [N II] /Ha for the whole galaxy is 0. 1 12±0.050. 

Due to its relatively low metallicity (compared with local spirals), the [N II] line is too weak 
to reach a S/N of 5 in individual OSIRIS pixels, but is clearly detected. We therefore integrate 
the spectrum within three annuli (see Figure [1] right panel). The physical lengths of the three 
outer radii of each annulus in the source plane are rl = 0.72 ± 0.1 kpc, r2 = 2.34 ± 0.2 kpc, 
r3 = 4.5 ± 0.4 kpc. The [N ll] line is robustly detected at S/N > 5 for the inner two annuli and is a 
3cr detection for the outer annulus which we give as an upper limit. The final spectra for the three 
annuli are presented in Figure |2] 
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Fig. 2. — Integrated spectra for the three annuli defined in Figure [Tl Vertical dashed lines indicate 
the zero-velocity positions of [N ll] and Ha. The y— axis is the observed intensity in units of 10^^^ 
ergs s^^ cm~^ A~^. The annuli radii rl, r2, r3 are chosen to ensure a S/N ([N ll]) to be >5 for the 
inner two radii. Blue curves are the Gaussian fits of the spectra. Note that the outer annulus r2 : r3 
is a 3 cr upper limit for [N II]. The declining ratio of [N II] /Ha from the central to the outer disk 
infers a decreasing metallicity with galactocentric distance. 

The metallicities converted from the N2 calibrator for the three annuli are 



154 ± 0.04 for r < rl 
12 + log(0/H) = I 8.38 ± 0.05 for rl < r < r2 

< 8.05 for r2 < r < r3 {3a upper limit) 



(2) 



The "global" metallicity measured from the integrated spectrum of all three regions is 8.36±0.04. 
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Note that the errors of line flux measureme nts are smaller than the intrinsic 0.18 dex scatter perti- 
nent to the N2 calibration (|Erb et al.ll2006|) . The s ystematic uncertainty of the N2 calibration with 
respect to other line diagnostics is about 0.04 dex (|Kewley & Ellisonll2008[) . 



A linear regression with censored data yields a metallicity gradient for the central 4.5 kpc 
of ^ ^° AR = —0.16 ± 0.02 dex kpc~^, or -0.15 dex kpc~^ if we treat the last data point as a 
measurement rather than an upper limit. The true gradient may be steeper than this value given the 
3a upper limit in the outer annulus. 



5. Metallicity Gradient Compared with Local and Higher Redshift 



In Figure |3] we compare the metallicity gradient of Spll4 9 with galaxies in the loca l and 
higher redshift universe, as well as the model predictions from IPrantzos & Boissieii (|2000l) . To 
avoid systematic effects caus ed by different calibratio n indicators, we use the oxygen abundance of 
H II regions for comparison. iKewley & Ellisonl(|2008|) showed that metallicities can be compared if 
the same metallicity calibration is applied consistently across all sample galaxie s and H ll regions. 
We co nvert all metallicities to the PP04 calibration using the calibrations in Kewley & Ellison 
(l2008l) . 



Prantzos & Boissieii (|2000|) modeled metallicity gradients of spiral galaxies at ages t = 3, 7.5, 
and 13.5 Gyr, for different velocities (Vc) and spin parameters (A). Since Spll49 is at an age 
of ~4.3 Gyr and has Vc=212±50 km s^\ we adopt model grids with 1^ = 220 km s~^, t = 3, 
13.5 Gyrs, and A = 0.03, 0.05, 0.07 for our comparison. The model metallicities have been cal- 
ibrated to the PP04 metallicity diagnostic. Spll49 shows a steeper gradient than predicted by 
Prantzos & Boissien (|2000h for Vc = 220 km s^. However, the steep gradient of Spl 149 is broadly 
consistent with the Vc = 150 km s^^ model (the dash-dotted model line in Figure |3]). 

Metallicity gradients in local galaxies are strongly correlated with Hubble type, b ar strength, 
and rnerging events; with more flattened gradients for barred galaxies and merging pairs (IKewley et al, 
2010l : lRupke etaPboiOah . Therefore it is significant that the metallicity gradient of Spl 149 is sub- 
stantially steeper than local late-type gala xies (see Figure |3]). F or the local galaxy gradients, we 
adopt the isolated spiral control sample of iRupke et al.l (|2010b|) comprised of 11 isolated mid-to- 
late type local galaxies with a mix of bar-strengths. The median slope for the local control sample 
is -0.041 ±0.009 dex kpc~^, with a standard deviation of 0.03 dex kpc^^. Th e median slope of gra- 
die nts from the Rupke e t al. (2010) s ample is consistent w ith the values of I zaritsky et al. ( 1994 ) 
and lvan Zee et al.l (Il998h samples. ThelZaritskv et al.l (|l994|) gradient sample composed of 39 local 
disk galaxies, and the Ivan Zee et al.l d 19981) sample added another 1 1 local spirals. Local early-type 
galaxies have even shallower gradients ([Henry & Wortheylll999|) . 
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A similar (and even steeper) gradient has been reported in a recen t work of Ijones et al.l (|2010a|) 
for a dispersion-dominated lensed system within only ~1 kpc at z=2. Cresci et al. ( 201o|) recently 
reported reversed metallicity gradients (i.e., the core is less enriched than the outer disk) for three 
un-lensed LBGs at 3. They explained the inverted gradients by inflow of pristine cold gas to the 
galactic center. Due to their low resolution (~ 4-6 kpc), the LBG selection, and potential presence 
of mergers and gal actic winds which may contaminate metallicity gradient, it is not clear how the 
Cresci et al.l (|2010|) results relate to the work presented here. 



In the right panel of Figure |3] we show that the metallicity gradient of Spll49 is still con- 
siderably steeper than local galaxies when expressed in units of dex/(i?/i?25), where i?25 is the 
5— band isophotal radius at a surface brightness of 25 mag arcsec"^. We estimate the i?25 for 
Spll49 by fitting the isophotal profile of the source-pla ne F8 14W-band imag e assuming a color 
difference of U — B = —0.1 mag for Sc type galaxies ([Fukugita et al.lll995h . The resulting -R25 



for Spl 149 is 12±2kpc, c omparable to the range of R25 in our local control sample ([Kewley et al. 



l3 



2010l : lRupkeetal.ll2010br) 



To analyze the effect of size and luminosity on the metallicity gradients, in Figure |4l we plot 
the metallicity gradients (in units of R25) v.s. the absolute fi— band magnitude Mb for local late- 
type galaxies and Spll49. Spl 149 is comparable to the brightest local late-type disks, while "the 
Clone Arc" is much more luminous (Mb = —22. 12). We see from Figure|4]that the scaled gradient 
of Spl 149 is steeper than in typical local galaxies. 



6. Discussion 



The redshift range of z=\ — 3 i s the perijjd when star-formation, mas s assembly and chem- 



ical enrichment activity p e aks (e.g.. 



Fan et al. 



200 ll: Dickinson et al.ll2003l : IChapman et al.l 12005 



Hopkins & Beacomll2006l : IConselice et al.ll2007h . Metallicity gradient studies of galaxies in this 
redshift range therefore give important insights into disk galaxy formation. The picture that the 
metallicity gradient flatt ens with cosmic time is qualitatively consistent with the "insid e-out" disk 
formation scenario (e.g., Prantzos & Boissier 2000; Bouwens et al. 1997 : Benson 20ld ). This sce- 
nario predicts that in the early stages of galaxy evolution, the inner galactic bulge undergoes a rapid 
collapse with vigorous star formation at the center, building a steep radial metallicity gradient from 
the core to the outer disk. Violent gas dynamics and accretion events marked these early stages of 
disk formation. Subsequent radial mixing and infall cause the gradient to flatten over the following 
Gyrs, resulting in a weaker gradient on average in late-type galaxies today. 



In this work, we have presented an OSIRIS IFU study of a grand-design lensed spiral galaxy 
at z = 1.49. The 22 x magnification allows us to resolve individual Hll regions at an intrinsic 
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Fig. 3. — Left: metallicity vs. galactocentric radius. Red lines are the measurements for Spl 149 at 
z=l .49 from this work. The gradient within the central 4.5 kpc is —0. 16 ±0.02 dex kpc^^. Vertical 
red dotted lines show the annulus used to average/sum the spectra. Purple dashe d lines show 



the typi cal gradients of local isolated late-type galaxies, using the control sample of IRupke et al. 



(|2010b|) . The orange dotted line represents the mean gradien t of local early-type gala xies, which 
is typically ~ 3 ti mes shallower than local late-type galaxies ( Henry & Worthey 19991) • Blue lines 
show the work of Ijones et al.l (|2010ah for a dispersion-dominated lensed galaxy the "Clone Arc" at 
z = 2.0 who report a even steeper gradient within a size of ~ 1 kpc. It is interesting to see that two 
z > 1 galaxies h ave considerably steeper gr adients than local galaxies. Black lines are the model 
predictions from IPrantzos & Boissien (l2000|) for rotational velocity V'c=220km s^ at age t=3Gyr 
(dashed curves) and t = 13.5 Gyr (solid curves), with spin parameter A=0.03,0.05,0.07. The black 
dash-dotted line is the model grid with velocity Vc=150km s~^ at t = 3 Gyr and A=0.03. The steep 
gradients of Spl 149 and "Clone Arc" are broadly consistent with the shape of the Vc=150km s^^ 
model grids. Right: The same as the left panel, except that the x— axis is expressed in the scaled 
radius R25. 



resolution of 170 parsec. We measure the radial distribution of chemical abundances in three 
annuli and find that the metallicity decreases from a central value of 12-i-log (0/H)=8.54 to less 
than 8.05 in the outer disk. The metallicity slope is —0.16 ± 0.02 dex kpc^^ for the inner 4.5 kpc, 
much steeper than that of late-type galaxies of comparable luminosities in the local universe, even 
when expressed in R25 units. Our results are broadly consistent with the "inside-out" model of 
disk formation. 

Metallicity gradient observations for a larger number of star- forming galaxies at different 
redshifts are required to build a well-defined sample to form a solid observational picture of the 



- 10- 



in 

CM 

or 



o 

cn 
O 



1 .0 
0.5 

0.0 

-0.5 

-1.0 

-1.5 

-2.0 
-2.5 



i 



• RupkelO 
O Zaritsky94 
A vonZee98 
■ Spl 149 




-17 



-18 



-19 



-20 



-21 



-22 



Fig. 4. — Metallicity gradient (dex/i?25) vs absolute B band magnitude Mb- The red square 
represents the measurement f or Spl 149 at z=lA9. Black filled circles indicate the locations of 
the local late- type galaxies of iRupke et al. JlOlOb I). Black empty circles are from the sample of 
Zaritsky et al.l (L1994), where the metallicity gra dients are measured for 39 local disk galaxies. 
Black empty triangles come from the sample of Ivan Zee et al.l (|l998h . where the metallicity gra- 
dients are provided for another 11 local spiral galaxies. Spl 149 has significantly steeper gradient 
than local galaxies of comparable luminosity, even when expressed in units of scaled radius. 



metallicity gradient evolution as a function of cosmic time. Before the next generation of tele- 
scopes, we are aiming to secure a few tens of gravitationally lensed systems as a pioneer sample 
to study the gradient evolution of galaxies at high-z. Future instruments such as the JWST, TMT 
and GMT will improve the current resolution limits by ~ 1 order of magnitude, revolutionizing 
chemical gradient evolution studies at high-2;. 
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